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TORQUE REQUIREMENTS FOR ORIENTATION OF A SOLAR 


BRAYTON SYSTEM IN EARTH ORBIT 
By Donald L. Sovers 
Levis Research Center 


SUMMARY 


Several torques that must he considered in the orientation control of a 
solar Brayton pover system are analyzed. The torques are divided into three 
groups: external, internal reaction, and torque necessary to compensate for 

variations in angular velocity. For each group of torques, general equations 
are derived and numerical values of the torques are calculated for a typical 
10-kilovatt solar Brayton system in a 300-nautical-mile ecliptic orhit. 


Of the external torques, gravitational and aerodynamic on the collector are 
the largest. They are sinusoidal and have maximums of 2.5x10“^ and 1.65x10’“^ 
foot-pound, respectively. Torque due to solar radiation on the collector is 
several orders of magnitude smaller. 


The internal reaction torque, vhich is caused hy accelerations of the 
turhomachinery, is the largest torque considered. A torque as large as 0.82 
foot-pound is due to acceleration of the turhoalternator during a typical alter- 
nator load transient. The internal torque, however, does not misorient the col- 
lector, hut rather rotates it about its axis, provided the turhomachinery axes 
are alined parallel to the collector axis of symmetry. If the turhomachinery 
axes are alined perpendicular to the collector axis, torques of 4.3x10"^ and 
2.0x10“^ foot-pound for the turhocompressor and turhoalternators, respectively, 
are required to keep the collector alined within l/4^. 

Torque necessary to compensate for variations in angular velocity is shown 
to he negligible. This torque compensates for angular acceleration with respect 
to the Sun hy applying an equal angular acceleration to the power system about 
its center of gravity. 


INTRODUCTION 

Studies are being made hy NASA to determine the feasibility of developing 
a solar Brayton system for space electric power generation. This system uses an 
all-gas thermodynamic cycle with the Sun as the energy source. A large parab- 
oloidal collector reflects the Sun's rays onto an absorber located at the focal 



point. Heat is transferred from the ahsorher to the working fluid (inert gas). 
The gas expands through a turhine which drives an alternator. The gas is then 
cooled, compressed, and returned to the ahsorher. Further discussions of 
Brayton cycle systems are given in references 1 and 2. 

A major requirement of this system is the proper alinement of the collec- 
tor with respect to the Sun. Studies have indicated that the collector must he 
alined within approximately l/4^ to insure effective operation of the power 
system. In order to determine the size and type of attitude control system re- 
quired to keep the solar collector alined, the magnitude and direction of the 
distiorhing torques must he known. This report presents an analysis of torques 
encountered hy a solar Brayton system in Earth orhit. 

The torques have heen divided into three groups: external torque, inter- 

nal reaction torque, and torque necessary to adjust for variations in angular 
velocity. External torque is caused hy external forces acting on the power sys- 
tem, such as gravity, atmospheric drag, or light pressure. Internal reaction 
torque is caused hy the acceleration of moving parts in the power system. The 
third torque is that torque necessary to change the angular velocity of the col- 
lector about its center of gravity. As the collector revolves about the Earth 
and Sun, its angular velocity with respect to the Sun is constantly changing. 

The collector, therefore, must he rotated about its center of gravity at the 
same rate of change to keep it in proper orientation with the Sun. 

Torques caused hy the connecting apparatus between the power system and the 
space vehicle have not heen considered. No consideration has heen given to the 
space vehicle to which the power system would he attached. 

General equations are derived for each of the torques, and then specific 
dimensions of a typical 10-kilowatt solar Brayton system are used to obtain val- 
ues of the torques. In this analysis, it is assumed that the solar collector 
and the power system are rigidly attached. For the numerical calculations the 
solar Brayton system has a 30-foot-diameter collector and a radiator 30 feet in 
diameter and 5 feet long. The system weighs 1500 pounds and is in a 
300-nautical-mile ecliptic orhit. 


SYMBOLS 

A area 

a semimajor axis of ellipse 

h semiminor axis of ellipse 

Cj) drag coefficient 

c distance from center of Earth to center of power system 
d half length of cylinder 

e eccentricity of ellipse 
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F force 

g gravitational constant at Earth’s surface 

h radius of cylinder 

I moment of inertia 

i, unit vectors along the rectangular coordinate axes 

K pressure constant 

k constant 

L distance from edge of radiator to origin of coordinate system 

I distance from origin of coordinate system to center of gravity 

M total mass 

m mass 

N normal to surface 

n distance from Earth to Sun 

P pressure from diffuse reflection of air 

p angular speed of turbine 

R radius of Earth 

r, Z cylindrical coordinates 

r^^ distance from elemental mass to center of Earth 

S solar radiation pressure constant at Earth’s orbit 

molecular speed ratio 
s surface 

T temperature 

ambient temperature 
surface temperature 
t time 

U velocity of Earth vith respect to power system 
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V velocity 

Vq velocity of power system about Earth 

^ot tangential velocity component of 

V velocity of power system about Sun 
s 

tangential velocity component of Vg 

V volume 

x^y^z set of rectangular coordinate axes 

a angular acceleration 

^ mis orient at ion angle 

X angle "between velocity vector and Z-axis 

6 density 

T) angle "between radiation vector and normal to surface 

p magnitude of position vector 

T torque 

cp angle of collector axis with line through center of Earth 

X angle defining orhit position 

^ angle "between velocity vector and normal to surface 

0, angular velocity of "body with respect to Earth 

CO angular velocity of body with respect to Sun 

Superscript: 

denotes vector 

EXTERNAL TORQUE 

There are several types of external torque. Of these, three types are 
dominant and will be considered in detail: gravitational torque, aerodynamic 

torque, and solar radiation torque. 

There are some other small external torques that should be mentioned but 
which are not calculated. These torques arise from the Earth* s magnetic field. 
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cosmic ray ’bom'bardment^ and the gravitational fields of other celestial hodies. 
The reaction of the power system to the Earth’s magnetic field depends on the 
electrical circuitry which has not yet heen specified. By proper design this 
torque can he made small. The latter two torques, cosmic ray homhardment and 
gravitational torque of other celestial hodies, are of the order of 10“^ times 
that due to solar radiation and gravitational torque from the Earth (ref. 3). 

Meteoroid homhardment, another external torque, is not calculated in this 
report. 


Gravitational Torque 

A hody in orhit about the Earth experiences a force due to gravity and 
centrifugal force. If the body’s energy is dissipated, these forces aline the 
axis with least moment of inertia with the center of the Earth. 

The torque acting on a hody about its center of gravity is 

T = A p X (1) 

*^v 

where p is the vector from the center of gravity of the hody to the elemental 
mass where the differential force dF acts. The force acting on an elemental 

mass due to gravity is g^R^yr^^dm, where R is the radius of the Earth, dm 

is the elemental mass, and r^ is the distance from the elemental mass to the 

center of the Earth. The centrifugal force acting on the elemental mass is 

p 

r^n dm, where a is the angular velocity of the hody around the Earth. The 
magnitude of the combined force acting on an elemental volume is 

dP = ^ (2) 

and is in the direction of -r^^ (that is, a vector from the elemental mass to 
the center of the Earth). 

A homogeneous right circular cylinder, as shown in figure 1, is considered 
as a model for the power system. The radius of th e cylinder is h and the length 
is 2d. By inspection of figure 1, the vector -r^ is 

-r^^ = (Z sin cp - r cos 6 cos cp)i - r sin 0 j - (c + Z cos cp + r cos G sin cp)k 

( 3 ) 

where r, 0, and Z are cylindrical coordinates of the elemental mass, c is 
the distance from the center of gravity to the center of the Earthy cp is the 
angle between the axis of the cylinder and the vector c, and i, j, and k 
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Figure 1. - Model of power system (right circular cylinder). 


are unit vectors along the x, and z axes, respectively. 

When equations (2) and (3) are used, the differential force can he -written 
in vector form as folloT^s: 

Z sin cp - r cos 6 cos cp)i 

- r sin 0 j - (c + Z cos cp + r cos 6 sin cp)kj (4) 

The vector p from the center of gravity of the cylinder to the elemental mass 
is 

p = (-Z sin cp + r cos 9 cos cp)i + r sin 0 j + (Z cos cp + r cos 9 sin cp)k 

(5) 

Substituting equations (4) and (5) into (l) yields 
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where the range of integration is 0<r<h, O<0<2jt, -d<Z<d. The dis- 
tance r^^ can he expressed as a fimction of c, r, Z, 0, and cp: 


V ? ? 2 

c + Z + r + 2cZ cos cp + 2cr sin cp cos 0 

The first integral in equation (7) is an odd function of 0j therefore, 
when integrating over 0 from -jt to n, the integral is zero. The second 
integral also is zeroj this can he seen hy integrating Z from -d to d and 
0 from 0 to 2 tc. Therefore, the total gravitational torque acting on the 
hody is given hy 
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( 9 ) 


The Integrals of equation (8) are evaluated in the appendix. The inte- 
gration yields the following expression for the gravitational torque: 

T = I Mg ^d^ - I h^j(sin 2cp)3 

where M is the total mass of the cylinder. 

It can he seen from equation (9) that the gravitational torque is zero if 
the moment of inertia is the same about all axes (d^ = s/i h^). Also it can he 
seen that the gravitational torque is periodic. The period is one-half the or- 
bit period with the angle cp directly related to the angle defining orbit posi- 
tion. The maximum magnitude of the torque occurs four times in each orbit. 

In order to estimate the magnitude of the torque for the power system as- 
sumed, the following dimensions are used in equation (9): h = 15 feet, 

d = 9 feet, and M = 1500 pounds. When these values are substituted into the 
torque equation, the gravitational torque about the center of gravity of the 
power system in a 300-nautical-mile ecliptic Earth orbit is 
-2. 45xl0"^( sin 2cp)j foot-pound. 


Aerodynamic Torque 

At altitudes up to 500 miles, air density is great enough to produce a sig- 
nificant force on a moving body. When the net force does not act through the 
center of gravity, a torque is produced. The aerodynamic torque consists of two 
parts - torque due to the incidence of the air molecules and secondly, torque 
due to the reflection of the air molecules. In this analysis the air molecules 
are assumed to be diffusely reflected. 

The torque on a surface is 


T 



p X dF 


( 10 ) 


where p is the ^ctor from the center of gravity to the point where the dif- 
ferential force dF acts. The component of the differential force due to the 
incidence of the air is 

^incidence = S I ^ ^ T I - cos r k) (ll) 

where Cp is the drag coefficient, 6 is the air density, U is the velocity of 
the air with respect to the body, \|r is the angle between the velocity vector 
and the normal to_th_e surfac_e, y is the angle between the velocity vector and 
the z-axis, and i, j, and k are unit vectors along the x, y, and z axes, 
respectively (fig. 2). For this analysis, Cp is assumed to be 2. 
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The component of the differential force due to the diffuse reflection of 
air molecules is given hy 




Reflection 


= P dA 


JL 

|W| 


(12) 


where P is pressure and N is the normal to the surface. Above 100 miles 
altitude, the air density is sufficiently low so that collisions between mole- 
cules after reflection can be neglected. The value of the pressure P is 
found in reference 4 to be 


P 


2 y T . 


cos 




where Sp is the molecular speed ratio, Tp is the surface temperature, and 

is the ambient temperature. Por a given orbit where the altitude is a con- 
stant, let P = K cos ^ where K is a constant and 


K 


= 2 y Ti Sr 


The total differential force then, adding equations (ll) and (12), is 



X 



Figure Z - Collector and radiator. 


cos ^1; dA 
(13) 

The power system has two components, 
the collector and the radiator, that contri- 
bute most of the surface area. Only these 
components, therefore, are considered in this 
analysis. 

Collector . - The collector is shown 
schematically in figure 2. It is a parab- 
oloid with axes x, y, and z fixed in the 
collector. The angle between a line joining 
the rim and the focal point and the z-axis 
is 60^. The radius of the collector is 
15 feet, and since the slope of the surface 
at the rim in the x, z-plane is the 

equation of the surface is 

+ y^ - 30y/3z = 0 

The normal to the surface is the gradient of 
the surface. Therefore, for this collector 
the normal to the surface is the vector 


= ^ 


ap 


5U^(-sin r i- cos r k)+ K, 


IT 

|N1. 
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N = xi + yj - 15 Vs k (14) 

where i, and k are unit vectors along the x, y, and z axes, respec- 
tively. 

If the angle x (13)) is the angle between the incident air molecules 

and the z-axis, parallel to the x, z-plane,_^then the unit vector representing the 
direction of the incident air is -(sin r)i - (cos y')k. When 0^ < T < 60^, the 
surface striking the air molecules is concave. This means that some of the re- 
flected molecules will strike the surface at other points. However, the con- 
cavity of the collector is small enough that this effect can he neglected. 

The cosine \lr is 


cos t = (xi+ yj - 15\/5 k)» (-sin r i - cos r 

Vx^ + + 675 

_ -X sin X + 15 Vs cos y 
Vx^ + + 675 


Expressed in cylindrical coordinates 

-r cos 9 sin y + 15 V5 cos y 


cos i|r = 




(15) 


+ 675 


An elemental area is 


OA = — - — + 675 dr d0 

15^3 

Substituting equations (14), (15), and (16) into equation (13) yields 


(16) 


[• 


dP = |6U^(-sin T i - cos Tk) + 


K 




(r cos 6i + r sin 0 J - 15 Vs k) 


+ 675 


'] 


(■ 


15^3 


cos 0 sin y + r cos yjdr dO 




(17) 


The vector p of equation (lO) is defined as the vector from the center of 
gravity of the power system to each elemental area. The center of gravity of 
the power system is assumed to lie along the axis of the collector and is lo- 
cated at z = Z. Therefore, 


p = r cos 0i + r sin 0 j + 


(— - 

\30->/3 / 


(18) 
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Substituting equations (17) and (18) into equation (lO) yields 


T = 


I 


Kr sin G 


“6U r cos Y sin G - " I 15 -v /5 - Z + — — = II 1 

Vr^ + 675 


( 


1SV3 - I + 


— )]i 

30 V3/j 




+ l6U“(r cos G cos x JC + I sin 

30 Vs 


in 


Kr cos G 




r + 675 


( 


isVs - i + 


30 V3/ 


+ ( 6U^r sin 6 sin x^k 




(- 


15 Vs 


cos G sin Y* + r cos y)d0 dr 


■where the range of integration is 0 < r < 15 and 0 < 0 < 2 tc. Integrating 0 
from 0 to 2jt (noticing that only terms containing cos^e alone or which are 
independent of 0 will not vanish) yields 


T = 


15x/3 


sin X 


ft 


6U^ cos x(soVs - 2r^) 


- K 




r"^ + 675 


( 


15-V3 - I + 


30 Vs 


:) 


&r 0 


and vhen integration is performed over r, the aerodynamic torque exerted on the 
collector is given hy 


T = 


Vs 


8U 


112.5 (Vs I - (' 


+ (442 Z - 12 721 

J-O 


) sin 7 ^ j 


(19) 


For a "body in a 300-nautical-raile circular orhit, U is 4.71 miles per sec- 
ond, the air density is about 1.5x10“^^ slugs per cubic foot (ref. 5), and K 
is 2.02x10"^ pounds per square foot. Using these values in equation (19) and 
assuming the center of gravity is at Z = 10 feet, yield the aerodynamic torque 


T = (1.86X10”^ sin 2x - 2.03x10"^ sin x)j ft-lb 


-4 


( 20 ) 
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It is noted that this expression is valid only for -60° < T < 60°. When 
60° < y < 90°, the air molecules are striking partly the front and hack of the 
collector, and when 90° < T < 120°, the air molecules strike part of the hack 
of the collector. It is difficult to get an exact solution for these ranges 
of y. When 120° ^ y ^ 240°, the air molecules are striking the entire hack 
of the collector, and the formula for the torque is the negative of the torque 
derived previously for -60° < y < 60° (eq. ( 20 )); y is therefore replaced 
hy y', which is measured from the -z-axis (as shown in fig. 2, p. 9). When 
y = 90°, exactly half of the collector is homharded hy the air molecules; there- 
fore, a calculation for this special case will give some insight into the case 
for 60° < y < 120 °. 

Since the air molecules are striking the hack of the collector only, the 
normal to the surface is 

N = -r cos ©T - r sin 0 j + 15 >/3 k 


The cosine ^ is 


cos \|f = 


_ N • (-i) _ r cos 9 


|K| 




r^ + 675 


Substituting these values for N and cos ^ into equation (13) yields the dif- 
ferential force 


dF = 


— f-SU^i + 

15-</3 |_ 


K 




+ 675 


(- 


r cos 0 i - r sin 0 j + 15 


Vs k) 


r cos 6 dr d 0 


( 21 ) 


Substituting equations (18) and (21) into equation (lO) yields 


T = 


0^ / (rjiji2j./i5V3+ . z) 

Vs jf j[Vr2 + 675 V “Vs /. 

* - 


Kr cos 6 


+ 675 


( 


isVs + 


30 Vs 




(a 


+ l&U^r sin 


r cos 0 dr d 0 


Jt 

The range of integration in this case is - 2 < 6<2 and 0 < r < 15. 
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Integrating over 0 from to + ^ yields 



If the center of gravity of the power system is at Z = 10 feet, the torque in 
equation (22) is 

T = 0.50X10"^ J ft-lb r = 90° (23) 

It can *be seen from equations (20) and (23) that the aerodynamic torque is 

cyclic. If the center of gravity is at Z = 10 feet^ the torque is positive 

for 0 < r < 180^ and negative for 180^ < y < 360^. The torque is zero for 

Y* = 0^ and 180^. The maximum occurs near y = 4-5^ and is 1.65x10"^ foot-pound. 

Radiator . - The second structure to consider for aerodynamic torque is the 
radiator. It is located as shown in figure 2 (p. 9). The model for the radi- 
ator is a circular hollow cylinder 30 feet in diameter, 5 feet long, and located 
L feet from the origin. 

The external surface is considered first. The air strikes the part 

• 5 ^ ^ S 5 when the range of y is 0 <tS 180°. Since the surface is con- 

vex, the reflecting air does not reimpinge on the surface. The total differen- 
tial force acting on an elemental area is given in equation (13). The normal 
to the surface is 

W = -r cos 0i - r sin gJ (24) 

The velocity of the air with respect to the hody is 

U = U(-sin vT - cos irk) 

The cosine of the angle between these two vectors is 

U • 1 

cos t - |U||U| 

cos ijr = cos 0 sin T (25) 

An elemental area is 

dA = r d0 dz (26) 
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SulDstituting equations (24), (25), and (26) into equation (13) yields 

dF = [sU^(-sinTi- cos y'k) + K(-cos 0T- sin ® T")^ <3.® <3-z (27) 


The position vector from the center of gravity is 

p = r COS0T+ r sin 0j+(z-'Z)k (28) 

Substituting equations (27) and (28) into equation (lO) yields 


= /{[- 


SU^r 


COS 


X sin 6 + (z - l)K sin sj i 


+ [^- ( z - l)5U^ sin X - {z - l)K cos 6 + rSU^ cos x cos ©J j 
+ rSU^ sin x sin ©k|(cos 9 sin T)r d© dz 

Jt ^ 

where the range of integration is L<z<L+5 and ~ 2 S © < "g* When the 
integral is evaluated, the torque exerted on the external surface is given hy 

T =||^5(L - 2) + ■^J(-2&U^ sin y - K) + 5 r&U^ cos y|r sin T J 

(0 < r < It) (29) 


A portion of the interior of the radiator is exposed to the air except at 
X = 0 and y = Jt/2. The interior surface will he analyzed in two parts: first 

where 0 < y < 80.5° and 99.5° < y < 180°, and secondly where 
80.5° < y < 99.5°. In the first range of y, the air will he assumed to strike 
half of the interior surface; that is, jr/2 < © < 3 jt/ 2. The error will he less 
than 20 percent with this assumption. For the second range of y, the air 
strikes a projected area in the form of an ellipse with a major axis 30 feet and 
minor axis of 30 cot y feet. Torque due to reflected air is not calculated 
because it is smaller than incident torque, and since the interior surface is 
concave, a large portion of the air will he reimpinging on the radiator surface. 
The differential force acting on an elemental area is 

dF = 6U^(-sin yi - cos yk)(-cos 0 sin y)r d© dz (30) 


Substituting equations 



(28) and (30 


sin © cos y i 
+ r sin 


) into equation (lO) gives 


+ 

r 


|^-( z - I )sin y + r 
sin © k\sU^(-cos 0 


cos © cos yjj 
sin y)r d© dz 


(31) 


14 



3 ■■ ■’ - 

Evaluting the integral in equation (31) gives 

T = -8U^r sin r{[io(L - 1) + 25jsin y + *1 rjt cos (32) 

-where 0 < y < 80.5^ and 99.5^ < T ^ 180^. 

The second range of y vill no"W he considered (80.5^ ^ T S 99*5°). For 
this range of y the range of variable z in integral (31) is changed to 
L + 5 + 2r cot y cos 0 < z < L + 5 for 80. 5^ < y < 90^, and 

L < z < L + 2r cos Q cot y for 90^ < T < 99.50. Evaluating integral (31) 

-with this nev range yields 

T = -8U^jtr^(L - I + 5)sin r cos r 3 80.5° < r < 90° (33) 

T = &U^«r^(L - Z)sin x cos T 0 90°^ xS 99.5° (34) 

If the radiator is located such that L = 10 feet, 2 = 10 feet, and 
r = 15 feet, the total torque on it in a 300-nautical-inile orhit is (by adding 

eq. (29) to eqs. (32), (33), and (34), respectively) 

T = (-6.98x10”'^ sin^r - 0.60>d.0“'^ sin ft-lb 

0 ^ r < 80.5° and 99.5° < r ^ 180° (35a) 

T = (-3.49x10“^ sin^r - 16.5x10"'^ sin x cos x - 0.60x10“'^ sin ft-lb 

80.5° < r < 90° (35b) 

T = ^-3.49x10“^ sin^r + 16.5x10“^ sin x cos y - 0.60x10“^ sin ft-lb 

90° < y < 99.5° (35c) 

From the symmetry of the radiator it can he seen that the following relations 
exist for the range of y, 180^ < T < 360^: 

T = (e. 98x10"^ sin^y - 0.60x10“^ sin y^3 ft-lb 

180° < y < 260.5° and 279.5° < y < 360° (35d) 

T = (3.49x10“'^ sin^y + 16.5x10"'^ sin y cos y - 0.60x10“^ sin y^j ft-lb 

260.5° < y < 2700 (35e) 

T = (3.49x10“^ sin^y + 16.5x10"'^ sin y cos y - O.OOxlO"*^ sin y)j ft-lb 

270° < y < 279.5° (35f) 
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This torque can he expressed as a function of orbit position 
tor is properly oriented vith respect to the Sun. The angle 
to X as follows: 


X if the collec- 
X is related 


X = 90O-r 0<X< 18QO 


X = 2700 - r 


"ISQO < X < 36QO for collector 
0 < X < 36QO for radiator 


The relation between aerodynamic torque and orbit position for the collec- 
tor and radiator is shown in figure 3. For comparison, the previously discussed 
gravitational torque is also shown as a function of orbit position X where 

X = 36QO - cp 


Solar Radiation Torque 

When a stream of radiation falls on a body, a pressure is exerted on that 
body. If the radiation makes an angle r\ with the normal to the surface, then 
the pressure on the surface is S cos^t), where S is the radiation pressure. 
The cos^T) term occurs because the components of both the force and the area 
on which the force acts must be considered. If the radiation is specularly re- 
flected, the pressure is doubled since the change in momentum of the radiation 
is doubled. 
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The solar radiation torque is calculated for the solar collector and for 
the ahsorher. On the collector the radiation pressure is 2S cos^t] assuming 
100-percent reflection. On the ahsorher, the pressure is S cos^-q assuming 
100-percent absorption. The pressure constant S includes light pressure, 
which is 9. 4 x 10~Q pound per square foot, and solar wind pressure, which is 
3. 5x10“^^ pound per square foot (ref. 6). Since solar wind pressure is negli- 
gible, the value of S used is 9. 4x10“^ pound per square foot. 

Collector . - Vhen a paraboloidal collector is not alined with the Sun, a 
larger amount of pressure will be exerted on the segments that are more nearly 
normal to the Sun's radiation. A torque is produced by this unequal pressure'. 

The torque on the collector about the center of gravity is 


( 10 ) 


where t, is the torque and p _i^s the position vector from the center of grav- 
ity to the differential force dF. The magnitude of the differential force is 

dP = 2S cos^Ti dA (36) 

where S is the radiation pressure constant, r| is the angle between the inci- 
dent ray and the normal to the surface, and dA is an elemental area. Expres- 
sions for cos T] and dA were given by equations (15) and (16), where y is 
replaced by the misorientation angle p. Substituting these equations into 
equation (36) gives 



15 Vs sin^B cos^9 - 2r^ sin p cos p cos 0 + 15 V^ r cos^P 

dF = 2S ^ ! ^ dr d0 

+ 675 

The direction of the differential force is normal to the surface; the differen- 
tial force vector being 


-I 


p X dF 


dF 



dF 


When the previous equation for W is used (eq. (14)), the differential force 
vector becomes 


dF 


^r cos 9 i + r sin 0 j - 15 V^ k } 


dF 




+ 675 


(37) 


Substituting equation (37) into equation (lO), along with the previous equation 
for "p (eq. (18)), yields the torque expression 
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When integrating over Q, the torque becomes 



When integrating over r, the torque equation is 

T = 2S^-1390 + 48.4 sin 2p j 

If the value of 8=9. 4x10"® pound per square foot is used, the torque is 

“ = (-2.61X10"'^+ 0.91x10"® z)sin 28 3 ft-lb (38) 

This is the solar radiation torque about an axis perpendicular to the z-plane 
and passing through the point (0,0, Z). It is noted that the torque expression 
is valid only for 0 < P < 30^, since for p > 30° part of the reflected rays 
strike the collector's surface again. Hovever, if the collector is properly 
alined throughout the orbit, P is alvays small (less than l/4^). 

Equation (38) indicates that solar radiation torque is small even for 
large mis orient at ion angles. For example, -with the center of gravity at 
Z = 10 feet and vith misalinement of 5°, the solar radiation torque is still 
only 3.0x10"^ foot-pound, which is several orders of magnitude less than gravi- 
tational and aerodynamic torque. 

Absorber . - When the collector is misorientated, the solar radiation on 
the absorber will cause a torque about the center of gravity of the power sys- 
tem. The torque is calculated for a flat plate located at the focal point of 
the collector and perpendicular to the axis of the collector, rather than for a 
hemispherical cavity as the absorber actually is. 

In order to determine the torque, it is assumed that during misorientation 
the image is displaced from the axis of the collector by the distance 
7.5V3 tan P feet on the focal plane. The image is assumed to keep its original 
size and original light intensity distribution. These assumptions will limit 
P to 1^ or less. With these conditions the torque is the product of the net 
force on the focal plane at zero misorientation and the distance away from the 
collector axis, which is 7.5\/3 tan p feet. 

To find the net force, the force from a solar ray reflected from a differ- 
ential area on the collector is multiplied by the cosine of the angle between 
the direction of the reflected ray and the normal to the focal plane. The 
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cosine of this angle can be written as a function of the cylindrical coordinates 
of the collector. For a ray striking the collector at distance r from the 
collector axis, the cosine is (675 - r^)/(675 + r^). The net force from solar 
rays striking the collector at r is 


S2irr 


675 - 
675 + 


dr 


The total net force from the collector is 


F 



S2nr 


675 - r^ 
675 + r^ 


dr 


Solving for F 


and using the value 
F = 


S = 9.4xlO“S 
4.82X10"^ Ih 


pound per square foot give 


The torque then is the product of F and 7.5>/3 tan P: 

T = 6.27X10"^ tan 6 "j (39) 

The torque vector has the direction of the positive y-axis when p is as shown 
in figure 2 (p. 9). 


The total torque exerted on the collector-absorber unit is the sum of the 
individual torques. If p is always small (less than 5^), sin p and tan p 
can be replaced by p. Therefore the sum of the torques, adding equations (38) 
and (39) (where I = 10 ft), becomes 

T = (-3.40X10"^ 6 + 6.27X10“^ 6^3 
= 2.87x10-4 pj 


INTEEWAL TORQUE 

Internal torque is due to acceleration of the turbomachinery. For the sys- 
tem considered there are two such machines, a turboalternator and turbocompres- 
sor. Torques must be exerted on these machines to change their angular veloc- 
ity. Equal and opposite torques are exerted on the power system. The torques 
are equal to the product of the moment of inertia of the rotating shaft and the 
angular acceleration of the shaft. The torques considered are only those which 
occur in transients from design point operation; torques in startup are not cal- 
culated. 

The turboalternator will be considered first. If there is a load change 
on the alternator, there will be a transient in rotating speed before the speed 
control can reestablish the speed at the design value. Typically, the speed may 
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vary as much as 1 percent from design in 1 second. The angular acceleration a, 
then^ would be 12.57 radians per second^ for a design speed of 12 000 rpm. A 
moment of inertia of 0,065 foot-pound second^ is representative of a turbo- 
alternator sized for the collector and radiator used in the discussion of ex- 
ternal torques. The torque^ therefore^ occuring during a typical alternator 
load change is 0.816 foot-pound. 

The turbocompressor has a greater speed range; however^ the changes in 
speed are slower. The torque will be calculated for a 10-percent change in 
speed from design in 60 seconds. For this unit^ a design speed of 38 500 rpm 
and a moment of inertia of 0.003 foot-pound second^ are representative. The 
turbocompressor torque^ therefore^ is 0.02 foot-pound. 

The torque needed to turn the alternator and compressor and to overcome 
bearing friction is experienced only inside the turbine packages. The vehicle 
does not have a torque exerted on it unless the angular momentum of the turbine 
packages changes. 

It should be noted that the internal torques are much larger than the ex- 
ternal torques. Of more significance is the change in angular momentum of the 
turbomachinery. The maximum change in angular momentum from steady- state oper- 
ation is 1.20 foot-pound second for the turbocompressor and 0.82 foot-pound 
second for the turboalternator. 

However^ if the turbomachinery axes are alined parallel to the axis of 
symmetry of the collector and radiator^ the internal torque only causes the 
power system to rotate about the z-axis and does not cause misorientation. If 
the turbomachinery axes are not alined parallel to the axis of the collector^ 
then internal torque does cause misorientation^ and therefore^ of prime inter- 
est^ is the torque required to keep the collector alined within l/4*^. 

Consider the turbomachinery axes to be alined perpendicular to the collec- 
tor axis; during a speed change in the turbocompressor^ the maximum change in 
angular velocity of the power system is 0.31x10"^ radian per second^ which is 
the quotient of the change in angular momentum (1,21 ft-lb sec) and the moment 
of inertia of the power system about the y-axis^ as shown in figure 1 (p. 6), 
with dimensions Z = 9 feet, h = 15 feet, and M = 1500 pounds. The torque 
then is the product of the moment of inertia and the angular acceleration a, 
which is that acceleration necessary to change the angular velocity of the power 
system by 0.31x10”^ radian per second without misorienting more than l/4C>, 

The misorientation angle (3 is 



The angular velocity p is 

P = at 

Eliminating time t gives 
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(3 


1 £; 
2 a 


Solving for angular acceleration a yields 

^ - 1 p! 

" 2 P 

The maximum misorientation angle alloved is 

B = 4 = 0.00436 rad 

4 

When (3 = l/^^^ P = 0.31x10“^ radian per second and then a = 1. 10x10 radian 
per second^. The torque is 


T = la = 4.3x10"^ ft-lh 

Therefore^ if a torque of 4. 3x10”^ foot-pound is applied to the pover system, 
the collector vill not misorient more than 1/4^ as a result of a change of speed 
of the turhocompressor. The corresponding torque for the turhoalternator is 
2.0x10“^ foot-pound. 


TORQUE COMPENSATING FOR VARIATIONS IN ANGULAR VELOCITY 

When the power system is in orhit about the Earth, its angular velocity 
with respect to the Sun is constantly changing. In order to aline the collec- 
tor with the Sun, the angular velocity of the power system about its center of 
gravity must be the same as the angular velocity of its center of gravity with 
respect to the Sun. Therefore, a torque must be applied to the power system to 
change its angular velocity and keep it alined with the Sun. 

There are two effects to consider, the changing angular velocity of the 

power system as it orbits the Earth and the changing angular velocity of the 

power system as it orbits the Earth in orbit about the Sun. The maximum angle 

that the collector misorients due to an orbit about the Earth can be shown to 

be small. If the collector is assumed to have a fixed orientation in space, 
the maximum misorientation angle p is 



where 2c is the diameter of the orbit about the Earth and n is the distance 
from the Earth to the Sun. For a 300 -nautical-mile ecliptic orbit, p is 
20 seconds of arc. 

For the solar Brayton power system this misorientation angle of 20 seconds 
is considerably less than the allowed misorientation angle and requires no cor- 
rection. However, because there may be systems that require greater orientation 


21 



Sun 



y 



accuracy than 20 seconds, and “because the variations in the angular velocity of 
the Earth orbiting the Sun cause the collector to znisorient more than is allow- 
able for the power system, a general discussion of these two effects of vari- 
ations in angular velocity follows. 

The velocity vector V of the power system in inertial space consists of 
two parts: the velocity vector of the power system about the Earth and 

the velocity vector Vg of the Earth about the Sun (as shown in fig. 4). The 
angular velocity of the power system with respect to the Sun o) is the tangen- 
tial component of the velocity divided by the distance from the Sun to the 
power system: 


Vot 

where n is the distance from the Earth to the Sun. The tangential velocity 
component of Vq is Vot = “^o Therefore, the angular velocity oa 

from equation (40) becomes 


-Vq cos X 
CD = + 

n 


' st 


n 


(41) 


where Vg-^/n is the angular velocity of the Earth with respect to the Sun. The 
angular velocity of the Earth with respect to the Sun (Vg-|^/n) can be derived 
from Kepler’s law, which states that the radius vector h sweeps out equal 
areas in equal time (see fig. 5 for definition of n and X^_). Therefore, 


~ dX^ = k dt 
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dX^ 

IF' 



(42) 


vhere k = jtab unit area per year. In this equation dX^/dt is the angular 
velocity of the Earth with respect to the Sun which is equal to 

In polar coordinates the equation of the Earth *s orbit as shown in fig- 
ure 5 is 


(n cos + ae)^ n^ sin^X^^ 

^ ^ — 1 


where e is the eccentricity of the ellipse. Solving for n gives 

2, 


ab 


n = 




e cos 


’‘i) 


b^ cos^X]_ + sin^X^^ 
Substituting equation (43) into equation (42) yields 


(43) 


dX^ 


27t 


(l - e2) 


3/2 


(l + e cos X^^)^ rad/yr 


(44) 


When equation (44) is substituted into equation (41), the total angular velocity 
with respect to the Sun is 


O) = 


-Vq cos X 


n 


— — ^ (l + e cos X^)2 


(45) 


(l - e2)^ 

If equation (45) is differentiated, the angular acceleration a is 


. V„(sin X)X sjt^e 


a = CD = 


(l - e2)^ 


2 

(^1 + e cos X^") sin X^ 


For a collector in a 300-nautical-mile ecliptic orbit, Vq = 4. 71 miles 
per second, X = 0.00109 radian per second, e = 0.0167, and n = 93 000 000 
miles. Therefore, the angular acceleration with respect to the Sun is 

a = 5.52xl0”ll sin X - 1. 32x10 ”1® ■*" ® ^1 ^s.6./ sec^ (46) 


The torque needed to turn the power system has two components: torque 

required to turn the mass of the system (t^_) and the gyroscopic torque required 
to turn the rotating turbomachinery masses (xg). The first component of torque 
is T-|_ = la, where a is the angular acceleration about the center of gravity 
which is equal to the angular acceleration of the power system with respect to 
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the Sim. The second term in the acceleration (eq. (46)) is negligible. The 
value of the moment of inertia I about the y-axis for the power system model 
shown in figure 1 (p. 6) is 3.9x10^ foot-poimd second^, if 1=9 feet, 
h = 15 feet, and M = 1500 pounds. Therefore, the torque is 

2.15x10”^ sin X foot-pound. 

The second (gyroscopic) component of torque is = Ipoo where I is the 
moment of inertia of the rotating mass and p is its angular speed. From 
equation (45), assuming e is zero, the angular velocity cd is 

0 ) = (-5.06 cos X + 19.92^10"® rad/sec 

The gyroscopic torque due to the turhoalternator is 

T 2 a = (-4.13 cos X+ 16. 27 ) 10 "® ft-lb 

and that due to the turbocompressor is 

T 2 c = (-6.11 cos X+ 24 . 1 ) 10 “'^ ft-lb 

If the turhoalternator and the turhocompressor are rotating in the same direc- 
tion^ then the total gyroscopic torque is 

Tg = (-4.74 cos X+ 18.68)10"® ft-lb 

The torques and ^2 perpendicular to each other. The torque 

is perpendicular to the ecliptic plane and Tg is in the ecliptic plane. 
Therefore, the torques must he added vectorially to get the total torque. When 
and Tg are added vectorially, it can he seen that the maximum value of 

the total torque required to compensate for variations in angular velocity is 
less than 0.01 of either of the large external torques (aerodynamic or gravita- 
tional) and is negligible. 

SUMMARY OF RESULTS 

Three types of torques exerted on a solar Brayton system in orhit were 
studied. They were external torque, internal reaction torque, and torque com- 
pensating for variations in angular velocity. 

Two external torques were found to he dominant, the gravitational and the 
aerodynamic torques. Both are periodic with gravitational torque having a max- 
imum of 2.5x10"^ foot-pound and aerodynamic torque having a maximimi of 
1.65x10”^ foot-pound for the collector and 0.75x10”^ foot-pound for the radi- 
ator. These torques will he at their maximum four times in each orhit. The 
solar radiation torque is several orders of magnitude less than the above 
torques even for large misorientation angles. 

The internal torques were the largest considered. The torque due to the 
acceleration of the turhoalternator may he as high as 0.82 foot-pound during a 
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typical alternator load change. The torque due to a typical acceleration of 
the turhocompressor may reach 0.02 foot-pound. These torques do not misorient 
the collector if the turhomachinery axes are alined parallel to the collector 
axis. If the turhomachinery axes are alined perpendicular to the collector 
axis^ torques of 4 . 3 xl 0“2 and 2.0x10-2 foot-pound for the turhocompressor and 
turhoalternator, respectively, are required to keep the collector alined 
within 1/4°. 

Torque compensating for variations in angular velocity vas found to he 
less than 0.01 of either of the dominant external torques and therefore is 
negligible. 


Levis Research Center, 

TTational Aeronautics and Space Administration, 
Cleveland, Ohio, August 27, 1965. 


25 



APEEEDIX - EVALUATION OF GRAVITATIONAL TORQUE INTEGRALS 


The first integral to evaluate from equation (8) is 

d y»h xi2tc 


/ 


/ 

-Zr d0 dr dZ 


1 

1 

+ Z^ + r^ + 2cZ cos cp + 2cr sin 9 cos 0) 


As 0 varies from 0 to 2jt the function varies from 


*2 

(c^ + + r^ + 2cZ cos cp + 2cr sin cp) 


(c^ + Z^ + r^ + 2cZ cos cp - 2cr sin cp) 

The average value of the function is about 



/ 2 2 2 

(c^ + + r^ + 2cZ cos cp j 

Therefore, integrating over 6, the integral becomes approximately 


// 


-Zr dr dZ 


+ Z^ + r^ + 2cZ cos cp) 


*'-d */0 

Integrating with respect to r yields 
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Since h is negligilDle compared to c, (c^ + + 2cZ cos cp^ 


becomes 


1/2 


/ 2 2 ^ 

\c^ + Z^ + 2cZ cos cp) = (c + Z cos cp) 


Then, 


2jt 


/■/ 


-Zr dr dZ 


(c^ + Z^ + r^ + 2cZ cos cp) 




-jth^ 


'-d ^0 

Integrating with respect to Z gives 
►d 


£ 


Z dZ 


(c + Z cos cp)^ 


-jth^ 




Z dZ 


ith^ 


(c + Z cos cp)" 


o 

cos cp 


[■ (c + 


cos cp) 2(c + Z cos cp)2 


d 

I 


Since d is negligible compared to c, (c^ - d^ cos^cp) becomes c^, then the 
first integral in equation (8) becomes 


2ith^ d^ 


cos cp 


The second integral to evaluate from equation (8) is 




r cos 0 d8 dr dZ 


/ 2 2 2 N 

Vc + Z + r + 2cZ cos cp + 2cr sin cp cos 0) 




2 2 2 
Since Z and r are negligible in comparison to c , the integral can be 

written as 


2 



1 

. \ 3/2 

2cZ cos cp + 2cr cos 0 sin cpy 


(c^ + 2cZ cos cp - 2cr cos 0 sin cp) ^ J 


r cos 0 d0 dr dZ 
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Expanding the tvo terms in the bracket by the binomial theorem and neglecting 
higher order terms enable the integral to be rewritten as follows: 


^d /*Jr/2 

t/-d Jo J 0 


-12 


3 ? 

r cos e 


-g sin cp d0 dr dZ 


Integrating with respect to 6 gives 


h /V/2 


-12 



r^ cos^e . j 

-T sin cp de dr dZ 


-5n . 

~ sin cp 


/'/ 


r-^ dr dZ 


Integrating with respect to r gives 


-3 -~r sin cp 
c^ 


t/-d Jo 


3 , 3jt h 

r dr dZ = — j j sin cp 


I 


dZ 


Integrating with respect to Z gives 


3jt ^ 

L 

c ■ 


- — — sin cp 


i: 


dZ = - ^ d sin cp 

2 ^4 ^ 


The torque in equation (8) becomes 

£ 

3 


T = 5gR^jt — d(s ^ h^^ sin cp cos cp j 


The mass of the cylinder is 


The torque^ therefore^ is 


M = Zrch'^dS 


T = -| Mg ^ (d^ - I h^)sin 2cp j 


( 9 ) 
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"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administratioyi 
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